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On) ■ Abstract 
> 

I We present the 0{a1) corrections to the cross section for the reaction 

CsJ I e~^e~ — > 7* — > QQ in the energy region close to the threshold. We assume 

that the energy of the reaction is such that both the perturbative expansion in 
Q\ • the strong coupling constant as and expansion in the velocity f3 of the heavy 

quarks can be used. We explicitly obtain terms 0(ag//3^, a^//3, a^) in the 
QhI relative correction to the threshold cross section. Using the ideas of asymp- 

CIhI totic expansions, we demonstrate how an expansion of Feynman diagrams in 

^ ' the threshold region is constructed. From this analysis we obtain a matching 



X 



relation between the vector current in full QCD and the quark-antiquark cur- 
rent in NRQCD at leading order in 1/m and the second order in the strong 
coupling constant. 

Theoretical predictions for the cross section of the reaction e+e^ —>■ QQ in the energy 
region close to the QQ threshold are of considerable interest for various phenomena. They 
are important for determinations of the b and c quark masses, as well as the value of the 
QCD coupling constant a;s(/i), with /i ~ 1 — 2 GeV, if one uses the sum rule approach for 
T and J/ip hadrons proposed in achievement of the accuracy is considered very 

important for these quantities. One should also mention the ongoing efforts to determine 
the decay rates of the heavy quarkonia to leptons with the 0{al) accuracy Also, for 
the future e^e~ or ^^fi~ colliders one is considering precision measurements of the top 
quark properties by studying its threshold production region. It is well known that for these 
and other applications, where the threshold region is of interest, the fixed order perturbative 
calculations break down and a resummation of the terms singular at threshold is mandatory. 
In the leading order, such resummation yields the well known Sommerfeld-Sakharov factor 
for the threshold production cross section. Aiming at the 0{al) accuracy for the threshold 
cross section, the first thing to be calculated is the perturbative expression for the cross 
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section at the order 0{a'^) in the energy region where -C /3 -C 1, where P denotes the 



quark velocity, (3 = \Jl — Am? / and s is the total energy squared 

From the technical point of view, the problem of perturbative calculations in the thresh- 
old region has never received much attention. On the one hand, it is clear that a small 
parameter in which a useful expansion could in principle be constructed — j3 — is there; 
on the other hand, it has not been quite clear how to use this parameter systematically in 
order to obtain a significant simplification of the Feynman integrals which should be cal- 
culated. On the conceptual level, the recognition of the existence of small and large scales 
in the threshold problems was formulated as the Non-Relativistic QCD (NRQCD), an ef- 
fective field theory which should be used to describe physics in the threshold region In 
this framework the corrections which originate from the scales ~ m can be incorporated 
through the so-called matching procedure. This requires calculations in the full QCD, in 
order to provide the matching conditions. It is expected that since the infrared behaviors 
of both QCD and NRQCD are the same, the matching calculations could be considerably 
simplified. To the best of our knowledge, however, it has not yet been demonstrated, how 
the matching calculations in QCD/NRQCD should be organized in order to achieve such 
simplification in practice. 

For practical reasons it is useful to be able to use the smallness of the relative velocity 
of the quarks on the level of individual Feynman diagrams, i.e. to formulate a prescription 
which operates with diagrams and subgraphs, rather than with composite operators and 
effective field theories. The advantages of such approach are its transparency and better 
control over the calculation. Indeed, similar approaches in other kinematic situations have 
recently permitted to complete many previously impossible calculations, both in QCD and 
in the electroweak theory (see for example 

In this paper we consider a variant of the asymptotic expansions which can be used in the 
threshold region. It allows to construct an expansion in /? of a given Feynman diagram. For 
the resulting Feynman integrals one can construct algorithms which permit their calculation 
in any order in /3, and can be encoded in a symbolic manipulation language. 

An approach of this type has recently been discussed in ||^. Although its correctness has 
not been proven, its construction is analogous to the well established asymptotic expansions 

Using this approach we calculate the 0{a1) correction to the cross section of the reaction 
e^e~ — i> QQ, up to terms of order relative to the Born cross section. In the result 

the terms with a color factor C|^, as well as the corrections induced by vacuum polarization 
insertions due to heavy and light fermions are identical to the previously obtained results in 



the Abelian gauge theory . This provides a non-trivial test of the asymptotic expansion 
method applied in the threshold region. The non-abelian terms presented below are a new 
result. 

We begin with introducing some notations. The cross section of the reaction e^e~ — *• 
7* — > QQ is written as 



TT / \ TT 



(1) 



where 
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s is the total energy squared of the reaction, m is the mass of the heavy quark Q, Nc = 3 
is the number of colors and eg is the charge of the quark Q in units of the electric charge. 
We also always use = as{m) in the MS scheme and the on-shell renormalization for the 
heavy quark propagators. 

The terms A*^^-* and A^^-* represent, respectively, the one- and two-loop corrections. 
Vertex corrections which contribute at these orders are shown in Fig. The term A*^^) 
is known in exact form as a function of (3. Near the threshold we are interested in the 
expressions for both A*-^'^^ up to and including terms We note, that up to this order 

in (3, there is no need to consider radiation of real gluons in the process of interest. With 
this restriction the expression for A*^^^ reads 



The two terms in this expression are known to be of rather different origins. The term 
proportional to the inverse power of f3 is the so-called Coulomb correction, while the second 
term is known to be a hard correction. Let us demonstrate how these corrections could be 
calculated using asymptotic expansions. 

In this approach for each diagram one has to identify those regions of the integration mo- 
menta which can give non- vanishing contributions at the given order in j3. In the calculations 
close to threshold, one should distinguish four different regions (the following description 
applies to the center of mass frame of the QQ pair) 0: 

1. hard region, where the momenta of the quanta are of the order of ~ m, 

2. soft region, where ko ~ |k| ~ m/3, 

3. potential region, where ko ~ and |k| ~ m/3, 

4. ultrasoft region k ~ 

One should construct all possible subgraphs of a given graph assigning the above labels to 
the lines in all possible combinations. After such assignment, the routing of the momenta 
should satisfy the "scale conservation." Namely, two ultrasoft lines, for instance, can not 
produce a potential line. On the contrary, two potential lines can produce an ultrasoft line 
and so on. 

To calculate A*^^^ to necessary order in /3, one has to consider the real part of the one-loop 
correction to the vertex 'y*QQ. In this diagram only potential and hard regions contribute. 
The other two regions generate massless tadpoles which vanish in dimensional regularization, 



^The term 0{P) in A^^) can also be used for the matching on non-relativistic resummed cross- 
section. This term represents a kinematical correction and leads, essentially, to the replacement 



^ = 7775 -4 + 



(3) 



7rV(2/?) 



■2(l + /?2)/(2/?), see 
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which we use throughout this calculation. The potential subgraph gives rise to a finite 
contribution 7r^/(2/5). The hard subgraph (which to the order 0{f3^) is obtained by simply 
considering the vertex correction at the point s = 4m^ in dimensional regularization) gives 
the constant term —4, if combined with the ultraviolet renormalization of the Born cross 
section. 

To see the importance of this classification we recall the following fact. The above result 
for A*^^) is valid in the region P ^ as, where perturbative calculations are still justified. 
When one approaches the region of small (3, one should perform a resummation of all terms 
of the form Such resummation results in the well known Sommerfeld-Sakharov 

factor, which is the modulus square of the fermion wave function at the origin, when the 
Coulomb interaction between nonrelativistic Q and Q is taken into account exactly. It is 
also believed that part of the subleading terms of the form a""'"^//?"' can be resummed by 
multiplying the Sommerfeld-Sakharov factor by the so-called hard correction: 



a = a(°)fl-4C^^)|^(0)|2 



It is interesting that this hard renormalization constant for the Sommerfeld-Sakharov factor, 
which has been known already for a long time, is the same as the contribution of the hard 
subgraph of the one-loop vertex correction and the ultraviolet renormalization. This is of 
course not accidental. If the approach based on asymptotic expansions in the threshold 
calculations is to be trusted, it is fairly clear, that the above formula is correct and the 
one-loop hard correction does indeed provide the renormalization of the Coulomb ladder to 
all orders in the coupling constant. A simplest way to see this is to say that the contribution 
of the hard subgraph provides a renormalization of the NRQCD vector current in order to 
match it on the complete vector current in QCD (see a more detailed discussion below). 

One can hope that the knowledge of the contributions coming from hard subgraphs at 
order O(a^) would permit a determination of the two-loop renormalization of the Coulomb 
ladder. 

We now turn to the consideration of the second order correction, described by the term 
A'^^\ It is convenient to decompose it into terms proportional to various SU(3) color factors 
(in SU(3) Cp = 4/3, Ca = 3, Tr = 1/2, and N^^h are the numbers of quark fiavors of mass 
and m, respectively): 

A(2) = Cf + Ca A^na + NlTr Af + NhTr Ag\ (5) 

(2) 

The only unknown term in this expression is A]yj^, which arises in the non-abelian theory. 

(2) (2) 

^A H L same in the abelian and non-abelian theory. In the framework of QED A)^ 

were obtained in |Q in the threshold region. The terms which describe the contribution of 

(2) 

both massless and massive fermions A^, jj were calculated in an analytical form for arbitrary 



f3 in |12|. Also, we would like to mention, that the leading 0{al/l3) terms in A^]^ were 



(2) 

predicted in |]13[ using the observation that similar terms in A^ are directly related to 



the contribution of light fermions to the Coulomb-like QCD interaction potential between 
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heavy quark and anti-quark. These terms can be also extracted from the expression for the 
cross section for e^e~ — QQ, which includes the resummation of the terms a^"'"'"^^/^"- 
calculation confirms these results and allows a determination of the terms in 

which is the main new result of the present paper. 
Various contributions to A*^^) are: 

4 2 4 

^ 12/52 (3 6 

+.^(_|_|,og,+ l.og2)+f-C3. (6) 

2/179 , ^ 8, \ 151 13^ _ 
+x^(^-logfl--log2)----C3. (7) 

' /3 VS 18/ 9 ^ ' 

Ag' = - - (9) 

^99 ^ ^ 

The terms i h coincide with the results obtained in QED. The non-abelian piece Aj^j^ 



was studied numerically in |jT3[|, where Pade approximation was used to obtain as a 
function of f3. A comparison of our result for A^^j^ with the results for the similar quantity 



presented in shows, that for /5 -C 1 there is a reasonable agreement. 



The details of the derivation of the above results cannot be presented here for lack of 
space. We only briefly explain how they are obtained using the asymptotic expansion in the 
threshold region. The contributions of all momenta regions, relevant for this calculation, are 
separated according to a classification given above and the integrands of the loop integrations 
are expanded in the respective small parameters (which vary from one momentum region 
to another). All divergences which arise in the course of this procedure are regulated using 
dimensional regularization. The calculation of the hard contribution is done by solving a 
system of recurrence relations which allows to reduce any hard integral to a limited set of 
master integrals. All integrals, where only potential lines are involved, are similar in their 
form to the integrals which one encounters in the non-relativistic perturbation theory and 
can be done easily. When some lines of a given diagram are either soft or ultrasoft, the 
calculation can be performed loop-by-loop. For this, the results for the one-loop eikonal 
integrals ITH] are useful as well as some results obtained in |TH . 



(2) (2) 

Let us elucidate the origin of various terms in the above expressions for and A^^j^. 
The most interesting contributions come from the hard subgraphs, because these contribu- 
tions could be used to discuss the renormalization of the Sommerfeld-Sakharov factor at the 
next-to-leading order as well as a matching of the NRQCD quark-antiquark vector current 
on the full QCD vector current at order Ola"^) and the leading order in 1/m. 

We note the following: in both AJ^ the terms which are not accompanied by powers 
of TT come from hard subgraphs. These are 39/4 — and —151/36 — 13(^/2, respectively, 
for A^'' and A^^j^. The terms of the order 0{l3^), which are multiplied by a single power of 
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TT^, are more tricky. First, one sees that there is a log (3 contribution in these pieces. This 
imphes that these terms get contributions from both hard and (all possible) soft scales, and 
these contributions are not finite separately. This in turn means that at the O(a^) order the 
matching coefficient of the NRQCD vector quark-antiquark current is not finite any longer. 
Removing this divergence by using the MS renormalization of the low energy effective field 
theory, we arrive at finite (but scale dependent) matching coefficients. 
We write 

^7V = (l - 2Cf^ + c,if,)Cp (^)'^ b'h^'Xk]^, (10) 

where the operator at the LHS of the above equation is the quark-antiquark current in 
the full QCD, while the RHS represents the quark-antiquark current in the effective field 
theory (with iph and Xh being two-component spinors), multiplied by a Wilson (matching) 
coefficient. We note that to the order 0{as, 1/m^) the matching condition for the quark- 
antiquark vector current has been obtained in |T^. 
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The knowledge of contributions coming from hard subgraphs allows us to obtain this 
matching coefficient directly: 

C2(Ai) = CpCA + CaCnA + NlTrCl + NhTrCh, (11) 

4-ic3. (12) 
~ T^'' ^^^^ 

cM = (14) 
c/^(/i) = -^7r2 + ^. (15) 

We note finally that within the present approach it is possible to calculate also the higher 
order terms in the expansion in /3. However, already for the terms of order in A*^^) one 
should incorporate the real gluon radiation. The expansion of the real radiation and the 
phase space in this case can be obtained as a slight generalization of a similar expansion. 



discussed in [|15| in the context of heavy quark decays. 

In conclusion, we have calculated the O(a^) correction to the heavy quark production 
cross section in e+e" annihilation in the threshold region, assuming <C /3 <^ 1. A method 
of systematic expansion of the Feynman diagrams near thresholds in powers and logarithms 
of the quark velocity f3 enabled for the first time an evaluation of the non-abelian terms. 
At the same time, comparison of the abelian terms with previously obtained results gives 
us confidence in the correctness of the method. The results of the present paper can be 
further used for matching-like calculations, to arrive at the predictions for the threshold 
cross section at the energy region, where /3 ~ a^. We have also given a matching relation 



6 



between full QCD and NRQCD effective currents to leading order in 1/m and second order 
in the strong coupling constant. 
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FIGURES 




(g) (h) (i) 



FIG. 1. Two-loop QCD corrections to the quark pair production; (a) is the one loop vertex 
correction; (b-i) are the two loop gluonic and fermionic corrections. 
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